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We study the superfluid pairing in a two-species gas of heteronuclear fermionic molecules with equal density. 
The interplay of the isotropic s-wave interaction and anisotropic long-range dipolar interaction reveals rich 
physics. We find that the single-particle momentum distribution has a characteristic ellipsoidal shape that can 
be reasonably represented by a deformation parameter a defined similarly to the normal phase. Interesting 
momentum-dependent features of the order parameter are identified. We calculate the critical temperatures of 
both the singlet and triplet superfluid, suggesting a possible pairing symmetry transition by tuning the s-wave or 
dipolar interaction strength. 
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I. INTRODUCTION 

The recent experimental realization and coherent control of 
high phase-space density quantum gas of the polar molecules 
40 K 87 Rb til-OA provides an excellent opportunity to study the 
effects of anisotropic long-range dipole-dipole interactions. 
Theoretical proposals employing degenerate polar molecules 
range from the study of exotic quantum phases of matter 1 4, 5 ] 
and quantum gas dynamics Ja, |7[] to quantum simulations of 
highly correlated condensed matter systems |8[] and schemes 
for quantum information processing [9]. 

Two fundamental properties of dipolar Fermi gases are su- 
perfluid pairing lflol4l2tl and Fermi surface deformation (TH- 
[Tal . originating from the partially attractive nature of the dipo- 
lar interaction and anisotropic Fock exchange interaction. For 
dipolar Fermi gases with two hyperfine states, one can tune 
not only the dipole-dipole interaction by a fast rotating ori- 
enting field flrjl . but also the s-wave inter-species interac- 
tion via a Feshbach resonance. Therefore, one expects that 
rich physics will emerge as a result of the interplay of the 
anisotropic long-range dipole interaction and short-range s- 
wave interaction. 

In this work, we study BCS pairing by taking account of 
the Fock exchange term in a self-consistent way. We find 
that the anisotropic nature of the dipolar interaction leads 
to an anisotropic momentum space distribution of the num- 
ber density and an anisotropic order parameter. We general- 
ize the definition of the deformation parameter introduced in 
Ref. II 1311 to describe the anisotropic number distribution in the 
pairing phase and find that it gives a good description. Inter- 
esting features of the order parameter in momentum space are 
revealed, manifesting fascinating consequences of the dipolar 
interaction. Competing effects of the contact s-wave interac- 
tion and the dipolar interaction are identified in the study of 
the transition temperature of the superfluid state, suggesting 
the possibility of tuning the pairing symmetry by tuning the 
dipolar interaction. 



ity, we further assume that each species has the same mass 
and density. The electric dipoles of the molecules with mo- 
ment d are oriented along the z— axis by a sufficiently strong 
external electric field such that the spin-independent part of 
the electronic dipole-dipole interaction becomes Vd<j(q) = 
(47r/3)d 2 (3 cos 2 6> q — 1), with 6> q being the angle between mo- 
mentum q and the direction of the z axis in which the dipoles 
are aligned. In addition, we assume that molecules also inter- 
act via a contact interaction with strength g. This system is 
described by the following Hamiltonian 

H — fin = ^(e fc - M)cL c k<r + 

kcr 

^ ^'(q) C k+q<T C P -q<T' C P CT ' Ck ^' (1) 

kpqcrcr' 

where /i is the chemical potential, n is the total number den- 
sity, V is the volume, and = fc 2 /2m (where we have set 
h = 1). The interaction potential V aa i (q) = <7<$<T.- ff '+Vd<i(q) 
contains both dipole-dipole and contact interactions. Antici- 
pating the importance of the Fock exchange term, we decou- 
ple the interaction in all three channels [17]: direct channel, 
exchange channel and Cooper channel, resulting in the fol- 
lowing effective mean-field Hamiltonian 



HmF — 2_^t tkcrC k(T Ckcr 
k 



kcr 

C Lr C - 



— kcr/ 



(2) 



Here £ kcr = £k — M + 9 n /2 + ^ka, with self-consistent mean- 
fields defined as 

£ kCT = --^V drf (p-k) <c pff c pCT >, (3) 
p 

A CT/ff (k) = - ^ V aa , (k - p) < c_ kCT ,c kCT > . (4) 



II. MODEL 

We consider a homogeneous gas of two species of 
fermionic heteronuclear molecules, a =t and J.. For simplic- 



Some comments are in order: The contact interaction affects 
the single-particle spectrum by shifting the chemical poten- 
tial, which may be redefined as fi = /i — gn/2. The self 
energy encodes the anisotropic dipolar contribution from 
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the Fock exchange term to the dressed single particle spec- 
trum, which justifies our treatment. In addition, both parts of 
the interaction contribute to the pairing field. 

The Hamiltonian © is diagonalized by invoking the Bo- 
goliubov transformation 111 811 . We obtain self-consistent equa- 
tions for the self energy, the order parameters and the number 
density 



]T^ d (p-k) 
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where E ka = + So-' |A ffff '(k)| 2 is the quasi-particle 

spectrum and /? = 1/fegT is the inverse temperature. Equa- 
tion (|6]l formally diverges. For the contact interaction, one 
can eliminate the interaction strength g in terms of the s-wave 
scattering length a s using 1/g = m/ (47ra s ) — A ^ k The 
dipolar interaction can be regularized by replacing the bare in- 
teraction V<ja' (k — p) with the vertex function T a a' (k — p) 
as explained in Refs. [ 10] and [ 19]. To first order in the Born 
approximation, the gap equations become 



A CTCT ,(k) = ^K CT '( k 
p 



tanh /3.EkCT 

2£ka 



1 
2lk~ 



A^(p). 



The above equation, together with Eqs. (O and (O comprise a 
complete description of the dipolar Fermi gas and needs to be 
solved self-consistently. 

Due to the symmetry of the interaction potential, the mo- 
mentum distribution, order parameter and self energy possess 
azimuthal symmetry in thermal equilibrium, as can be seen 
from the self-consistent equations by integrating out the az- 
imuthal degree of freedom tfik- Thus, the physical quantities 
rik, A crcr '(k) and Sko- are only functions of (fc,#k)- Numer- 
ically, we parameterize these quantities by two-dimensional 
grids living on the domain fi = [0, k c ] X [0, 7r], where k c is the 
momentum cutoff. For spin singlet pairing the order param- 
eter possesses inversion symmetry as well as azimuthal sym- 
metry, A(fc,6*k) — A(k,n — 6k), while for spin triplet pairing 
A(fc,6>k) = — A(fc, 7T — #k)- In our calculation, we parameter- 
ize the dipolar interaction by the dimensionless coupling pa- 
rameter Cdd = md 2 {n (T ) 1 / 3 ' ', where n a = kp/(6n 2 ) = n/2, 
and the Fermi energy is Ep — k F /2m. In addition to the az- 
imuthally symmetric solution, vortex states with azimuthally 
varying phase (arg A = ivip, where ip is the azimuthal angle 
and v is integer) are expected to exist and can be treated by 
straightforward generalization of the presented approach. 



III. SPIN SINGLET PAIRING 

To investigate the interplay of the contact and dipolar in- 
teractions, we devote in this paragraph to studying spin sin- 
glet pairing at zero temperature. In the normal phase, the 
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FIG. 1. (Color online) Deformation parameter a as a function of s- 
wave coupling strength a at Cdd = 1 for the singlet superfluid. 



Fermi surface of the dipolar gas has an ellipsoid shape [13]: 
n(k) = Q(k F — o?k 2 z — k 2 /a — k 2 /a), where a is the de- 
formation parameter. In the superfluid phase, we propose to 
use the deformation parameter a to measure the anisotropy 
of the single-particle momentum distribution with a simi- 
lar strategy. The angular distribution ng(9k) can be ob- 
tained by integrating out the magnitude and azimuthal an- 
gle of the momentum ng(9k) = J n(k)k 2 dk/(2ir 2 ). Since 
ng(0)/ng(ir/2) = k\jk\ and a 2 k 2 = k 2 /a, we can de- 
duce that a = [n g (n/2)/ng(0)} 2/9 . In Fig. Q] the deforma- 
tion parameter a is plotted as a function of the s-wave cou- 
pling strength, characterized by l/kpa s , at fixed dipole inter- 
action strength Cdd = 1- As the s-wave coupling strength 
is increased, the deformation parameter a increases before it 
reaches a maximum, signaling that the anisotropy of the sin- 
gle particle momentum distribution is mitigated. It is energet- 
ically favorable for the system to become less anisotropic to 
benefit the energy gain from BCS pairing. The appearance of 
a pronounced peak suggests that there is an optimal value of 
the s-wave interaction strength to mitigate anisotropy. 

The single-particle momentum distribution is shown in 
Fig. |2] The angular distribution ng is symmetrical with re- 
spect to the orbital axis, with ng(9k) — ng(ir — 9k) ex- 
tending to the polar axis with a parabolic shape, as sug- 
gested in panel (a). The angle-averaged momentum distri- 
bution rik(k) = J n(\z)<Kl/ '4ir is shown in panel (b). For 
k < kp, the quantum states are essentially fully occupied. 
For k > kp the occupation number decreases rapidly with a 
fc~ 4 tail, as can be derived from Eq. (0. This fc~ 4 asymptote 
is valid for wave numbers smaller than the inverse range of the 
two-particle interaction. It is characteristic of BCS-type pair- 
ing and is also predicted in isotropic Fermi gases li20L 12111 . To 
further reveal the anisotropic nature of the momentum distri- 
bution we plot the radial momentum distribution for different 
polar angles 9k = 0, 7r/6, 7r/3, tt/2 in Fig. [3] For any an- 
gle 9k, the momentum distribution drops very sharply from 
almost full occupation to very small values at a critical value 
k c . At 9k = the value of k c is the smallest. As one in- 
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FIG. 2. (Color online) (a) Momentum integrated angular num- FIG 4 (Color online) Radial gap parameter profile in mome ntum 

ber distribution; (b) Angle-averaged momentum number distribution space at 1/kpa = _ x and Cm = 1 and different polar angles for 

nk (k) = J n(k)dfi/47T at 1/fe.Fa = -1 and Cd d = 1 for the singlet me smg i et SU p e rfluid. 
superfluid. 
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FIG. 3. (Color online) Radial number distribution in momentum 
space at 1/kFa = —1 and Cdd = 1 and different polar angles for 
the singlet superfluid. 



creases the polar angle from to tt/2, the critical momentum 
increases monotonically. The critical momentum fc c ((9k) as- 
sumes a similar role as the Fermi surface for the normal Fermi 
gas and defines an ellipsoidal shape with the long axis point- 
ing along the z axis and the short axis in the x — y plane. 
Thus, our definition of the deformation parameter a above is 
justified. 

Let us turn the attention to the anisotropy of the order pa- 
rameter, which reveals the excitation spectrum of the system. 
As shown in Fig.@] the order parameter shows interesting fea- 
tures in the momentum space: the order parameter is a non- 
monotonic function of the magnitude of momentum along 
fixed polar angle, a feature similar to the one found 111 111 in 
the spatial distribution of order parameter for trapped polar 
molecules. For different polar angles the magnitude of order 
parameter coincides at two typical momentum magnitudes, 



approximately and kp. Around the Fermi surface (k « ftp) 
the slope dA/dfc increases monotonically from negative val- 
ues for #k = to positive values for 0-^ = tt/2. The isotropy 
of the order parameter at the Fermi surface and at zero mo- 
mentum is a striking feature of the solutions of the mean-field 
Hamiltonian (2). 



IV. CRITICAL TEMPERATURE 

At T > T c , the gap parameter vanishes and the system un- 
dergoes a 2nd order phase transition to the normal state. For 
spin triplet pairing the order parameters A aa > (k) are decou- 
pled, satisfying the gap equations with exactly the same form. 
For convenience, to determine the critical temperature, we de- 
note the order parameter A aa r (k) by Ak, regardless of the 
pairing spin symmetry. As T approaches T c from below, the 
gap parameter is a small parameter, and the gap equation can 
be linearized: 



l/2e p 



tanh(/3|£ p |/2) 



We define 



ft(p,k) = V(p-k) 



l/2e p - 



tanh(/3|£ p |/2) 



then one gets an eigenvalue equation 



A k = ^ft(p,k)A F 



(8) 



Let us define the orthonormal function <fii(cos8) = 

^Y^-Pi(cos9) and expand the gap parameter in terms of 

this basis functions Ak = J2i A;(fc)</>/((9k). We can rewrite 
the Eq. © as 



/ p 2 dp 



5^fl«,pt'A,,(p), 
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FIG. 5. (Color online) (a) Critical temperature as a function 
of dipole-dipole interaction strength Cdd for singlet superfluid at 
1/kFO, = —1 and for triplet superfluid; (b) critical temperature as 
a function of 1/kpa at Cdd = 1 for singlet superfluid. 



where 



We show the calculated critical temperature in Fig. [3] In panel 
(a), for fixed s-wave coupling with 1/kpa = — 1, the critical 
temperature T c of the singlet superfluid decreases monotoni- 
cally as the dipole interaction strength increases. This is due 
to competing effects of the s-wave and dipolar interactions. 
While the s-wave interaction favors a spherical Fermi surface 
over an ellipsoid Fermi surface favored by the dipolar inter- 
action. In the case of spin triplet pairing, T c increases with 
the dipolar interaction as expected since the s-wave interac- 
tion does not participate in the pairing mechanism. When the 
system is cooled down from the normal phase, the actual pair- 
ing symmetry, spin singlet or spin triplet, is determined by 
which of the critical temperatures are higher. By tuning the 
dipolar interaction strength Cdd, the quantum phase transition 
between singlet and triplet superfluids may be realized. In 
panel (b), for Cdd = 1, the critical temperature for the sin- 
glet superfluid is clearly a monotonic function of the s-wave 
coupling strength. Here we only show T c in the range of weak 
interaction where the mean-field description gives quantitative 
reasonable results. For strong interaction at finite temperature, 
the fluctuation contribution is significant and one needs to re- 
sort to methods beyond the mean-field description. 



J d(cos6*k) J d(cos# p )7£(p — k)0/(cos0k)<^'( cos ^p)- 

(9) 

In matrix notation we may write 

A = KA, (10) 

p 2 Sp 

T^kLpl' — /n s 9 Ttkl,pl>, (11) 

( 27r ) 

£ w =A,(ft). (12) 

In this equation, the matrix and vector indices are jointly la- 
beled by the momentum magnitude k and the orbital index I. 
For singlet pairing I runs over even and positive integers while 
for triplet pairing I is odd. The critical temperature T c can be 
found from the solutions of the secular equation 

Det(£-I) = 0. (13) 



V. EXPERIMENTAL SIGNATURE 

The anisotropic nature of the superfluid order parameter 
and momentum distribution bear consequences for experi- 
mental observations. The single-particle momentum distri- 
bution of trapped Fermi gases is routinely observed by time- 
of flight measurements II22I1 . Collecting angle-resolved data 
from a gas released from an axially symmetric trap should 
reveal the anisotropies predicted in Section III. Anisotropic 
features of the pairing gap and further details of the single- 
particle s pec trum could be probed by radio frequency spec- 
troscopy miMi- 
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